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Coordination sequences of periodic and quasiperiodic graphs are analysed. These 
count the number of points that can be reached from a given point of the graph 
by a number of steps along its bonds, thus generalising the familiar coordination 
number which is just the first member of this series. A possible application to the 
theory of critical phenomena in lattice models is outlined. 

1 Introduction 

The locations of critical points of lattice models (such as Ising models, per- 
colation problems or self-avoiding walks) depend on the topological structure 
of the underlying graph, in particular on the dimension and on the (mean) 
coordination number. 1-5 However, this can only be a first approximation, and 
more detailed information on the graph is required for a better understand- 
ing of the precise location of critical points. For this purpose, quasiperiodic 
graphs prove helpful because they provide different cases with the same mean 
coordination number and also a check of the "universality" of this approach. 

Below, we display first results on the combinatorial issue to calculate the 
coordination sequences, and on general tendencies visible. To this end, we 
summarise the results for root graphs 6,7 and present the averaged coordination 
numbers for the rhombic Penrose and for the Ammann-Bccnker tiling. 

2 Coordination sequences 

A natural family of periodic graphs is given by the so-called root lattices and 
their relatives, the root graphs, 7 encoded by the Dynkin diagrams of A n , B n 
(n > 2), C n (n > 3), D n (n > 4), G2, Fa, E 6 , E 7 , and Eg. In all these cases, 
the coordination sequences s(k), where s(k) is the number of fc-th neighbours 
and s(0) := 1, can be calculated completely and encapsulated in an ordinary 



Table 1: First coordination numbers of the Penrose tiling [r = (1 + %/5)/2]. 



k 


s(k) 


num. value 


k 


s(k) 


num. value 


1 


4 


4.000 000 


6 


980 - 588r 


28.596 015 


2 


58 - 30r 


9.458 980 


7 


-1614+ 1018r 


33.158 601 


3 


-128 + 88t 


14.386 991 


8 


2688 - 1638r 


37.660 326 


4 


288 - 166t 


19.406 358 


9 


-3840 + 2400r 


43.281573 


5 


-374 + 246t 


24.036 361 


10 


4246 - 2594r 


48.819 833 



generating function of the form 

- = (i^r (i) 

where n is the dimension of the lattice (resp. the graph) and P(x) is an integral 
polynomial of degree n. For the (hyper-)cubic lattice Z™, this polynomial is 
simply P n (x) = (1 + x) n . For the other cases, the result is 

n 2 n 

p *» = E CP' p *» = E 

fe=0 fe=0 
n n 

^„(-)=E(2fc)- fe ' ^nW=E 
fe=0 fe=0 

the expressions for the exceptional cases can be found in Refs. 6 and 7. Many 
other cases can be calculated, in particular periodic graphs with more compli- 
cated fundamental domains, where an averaging procedure over the possible 
starting points is necessary. Details on this will be given elsewhere. 

In view of the applications, one is interested in a variety of graphs that 
share the same mean coordination number, but have different (averaged) higher 
coordination numbers. Such examples are provided by planar rhombic tilings 
with TV-fold rotational symmetry As they are solely built from rhombi, the av- 
erage coordination number is always exactly four, but the second (and higher) 
numbers may vary. To expand on this, we present the coordination numbers 
for two examples, the rhombic Penrose and the Ammann-Bcenker tiling (also 
known as the standard octagonal tiling). The values for the Penrose tiling 
given in Table 1 were obtained by an exact window algorithm, whereas the 
results for the Ammann-Bcenker case in Table 2 are conjectures based on the 
averaged coordination numbers of a periodic approximant with 275807 vertices. 
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Table 2: First coordination numbers of the octagonal tiling (A = 1 + v2). 



k 


s(k) 


num. value 


k 


s(k) 


num. value 


1 


4 




4.000 000 


11 


360 - 128A 


50.980 664 


2 


48 - 


16A 


9.372 583 


12 


828 - 320A 


55.451660 


3 


-24 + 


16A 


14.627 417 


13 


-508 + 236A 


61.754401 


4 


28 - 


4A 


18.343 146 


14 


-996 + 440 A 


66.253 967 


5 


52 - 


12A 


23.029 437 


15 


2580 - 1040A 


69.217 895 


6 


48 - 


8A 


28.686 282 


16 


1620 - 640A 


74.903 320 


7 


-324 + 


148A 


33.303 607 


17 


-5288 + 2224A 


81.210 963 


8 


732 - 


288A 


36.706 494 


18 


2372 - 948A 


83.325 543 


9 


380 - 


140A 


42.010 101 


19 


1324 - 512A 


87.922 656 


10 


-1140 + 


492A 


47.793 073 


20 


1196 - 456A 


95.118 616 



3 Critical points 

So far, only rudimentary results exist. Critical points of periodic systems have 
been investigated in some detail, but the results are not fully convincing, and 
they are still rather incomplete. Furthermore, most results concern hypercubic 
lattices and hence graphs with different dimensions and coordination numbers, 
which obscures the dependence on higher order coordination numbers. 

For the problem of self-avoiding walks, there exist results also on a variety 
of quasiperiodic tilings, 2 and their analysis supports our claim that higher order 
coordination numbers move the location. In particular, it is clearly seen that 
the critical point increases with growing second coordination number, while 
the critical exponents seem to be universal, as expected. However, in order to 
arrive at a conclusive statement, more, and in particular more precise, data on 
graphs with coinciding dimension and mean coordination number are needed. 

4 Outlook 

It is an open question whether one can find explicit expressions for the gen- 
erating functions S(x) for quasiperiodic graphs. Clearly, these will not have 
the form of Eq. (J|) — an analysis of the generating functions of periodic ap- 
proximants shows that the polynomial degrees of numerator and denominator 
grow with the size of the approximant. Still, the coordination sequences show 
quite a regular behaviour, compare Fig. 1, also when compared to slightly more 
complicated periodic graphs as, for instance, the Kagome and the diced lat- 
tice. The latter show some similarity to the (dual) Penrose tiling with regard 
to high-temperature expansions. 8 Their mean coordination numbers s(k) grow 
linearly with fc, but with two different slopes for even and for odd values of k. 
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Figure 1: Coordination sequences of some periodic and quasiperiodic graphs. 



Ideally, one would hope for a kind of integral transform from the generating 
function of the combinatorial problem to the partition sum of the corresponding 
model of statistical mechanics. Such a transformation, however, cannot exist 
in general since these models are usually not solvable, so their partition sums 
have an analytic behaviour incompatible with that of our generating functions. 

Therefore, one can only expect an empirical formula for the correction to 
the location of critical points due to higher order coordination numbers. To 
get a better understanding of this, we intend to analyse the critical structure 
of various models and on a variety of graphs in the near future. 
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